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THE SADDLE-NODE TRANSCRITICAL BIFURCATION IN A 
POPULATION MODEL WITH CONSTANT RATE HARVESTING * 

K. V. I. SAPUTRA 1 ", L. VAN VEEN*, AND G. R. W. QUISPEL+ 

Abstract. Wc study the interaction of saddle- node and transcritical bifurcations in a Lotka- 
Volterra model with a constant term representing harvesting or migration. Because some of the 
equilibria of the model lie on an invariant coordinate axis, both the saddle-node and the transcritical 
bifurcations are of codimension one. Their interaction can be associated with either a single or a 
double zero eigenvalue. We show that in the former case, the local bifurcation diagram is given by 
a nonversal unfolding of the cusp bifurcation whereas in the latter case it is a nonversal unfolding 
of a degenerate Bogdanov-Takens bifurcation. We present a simple model for each of the two cases 
to illustrate the possible unfoldings. We analyse the consequences of the generic phase portraits for 
the Lotka-Volterra system. 
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1. Introduction. Interactions between bifurcations of equilibria and of cycles 
occur naturally in dynamical systems with parameters. Often the interaction points 
act as organising centres in the bifurcation diagram. At such points curves of local 
and global bifurcations converge and the behaviour of the system is determined to a 
large extent. Consequently, these interaction points have been the subject of intensive 
research over the last decades and all interactions which occur generically in systems 
without special structure have been classified and parsed in the literature. 

The bifurcation theory for systems with a special structure is, as yet, incomplete. 
In systems with a special structure bifurcations can have a lower codimension than 
that in the general case. For instance, the presence of a Z 2 symmetry in the dynam- 
ical system can render the pitchfork bifurcation codimension-one. This happens in 
particular in certain normal forms with S 1 symmetry, such as the saddle-node-Hopf 
normal form, after decoupling of the angular variable \7\ . For this reason interactions 
with the pitchfork bifurcations have been extensively investigated. To mention a few 
contributions, Scheurle & Marsden [T3] particularly discussed the existence of tori 
and quasi periodic flows resulting from saddle-node-Hopf bifurcations, while Broer & 
Vegter [4] discussed the existence of Shilnikov bifurcations. The existence of hetero- 
clinic orbits was investigated by Lamb et ai. [12] for the saddle-node-Hopf system 
with time reversal symmetry and for the saddle-node-pitchfork system by Kirk & 
Knobloch [10]. 

Interactions with the transcritical bifurcation, in contrast, have not attracted 
much attention. A reason might be that this bifurcation is not associated with a 
global phase space symmetry in contrast to the pitchfork bifurcation. Transcritical 
bifurcations, however, appear frequently in applications, for instance in predator- 
prey interactions [5], in mathematical models for the spread of diseases [5] or as a 
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model for phase transitions in plasma physics [2J. For the analysis of these models, 
it is important to know the dynamics organised by interactions of the transcritical 
bifurcation with other local bifurcations. In the current paper we will investigate the 
interaction between saddle-node and transcritical bifurcations. To our best knowledge, 
this interaction has never been reported on in the literature before. 

One simple setting in which we can see a codimension one transcritical bifurcation 
is that of a planar system which, possibly after a change of coordinates, has an invari- 
ant manifold which coincides with a coordinate axis, independent of the parameters. 
An example of a class of models with this property comes from population dynamics. 
In Lotka-Volterra type models the variables are the population densities of several 
species. If a species dies out it cannot be regenerated and therefore the coordinate 
axes in such a model are invariant and the origin is always an equilibrium state (see, 
e.g. [T7]). Here, we will focus on a Lotka-Volterra model that has been modified to 
include a constant term, which represents harvesting or migration. 

From a bifurcation theory point of view the transcritical bifurcation can be con- 
sidered as a nonversal unfolding of the well-know saddle-node bifurcation. The saddle- 
node bifurcation has the normal form 

x = fi + x 2 , 

whereas the normal form of the transcritical bifurcation is given by 

y = ay + y 2 . 

If we apply the transformation z = y + ~ to the system above we obtain 

« 2 2 
Z = -~ +Z > 

which is a normal form of saddle-node bifurcation parametrised by a. Thus, we can 
consider the transcritical bifurcation as an unfolding of the saddle-node bifurcation. 
Because the map /i = — a 2 /4 is non invertible at the bifurcation point a = fi = 0, this 
unfolding is nonversal. 

Using the idea above, we investigate two different SNT interactions, corresponding 
to a single and a double zero eigenvalue. In the former case, no additional bifurcations 
take place and the bifurcation diagram around the interaction can be obtained as a 
nonversal unfolding of the cusp bifurcation. The second case is more involved. The 
normal form of an equilibrium with two zero eigenvalues is the Bogdanov-Takens (BT) 
normal form. However, due to nondegeneracy conditions of the transcritical bifurca- 
tion, we obtain the normal form of degenerate BT (DBT) bifurcation. In addition 
to the saddle-node and transcritical bifurcations, Hopf, homoclinic and heteroclinic 
bifurcations appear. We find two topologically different diagrams corresponding to 
different unfoldings of the DBT singularity, named the elliptic and the saddle case [6] . 

We illustrate all different SNT interactions with a Lotka-Volterra model with 
a constant term, which can be thought of as constant rate harvesting or migration 

mm- 

±i = xi(bi + anxi + ai 2 x 2 ) + e, 

X2 — x 2 (b2 + Qa\X\ + a 2 2X2)- (1-1) 

Without the constant term, the origin is an equilibrium and both the cci-axis and the 
X2-axis are invariant. With the constant term included this equilibrium is displaced 
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Fig. 1.1. Bifurcation diagram of il.lV with e and 02 as parameters. The two saddle-node- 
transcritical interactions have been marked STi and ST2. Top: the saddle case for ST2, with 
b\ = 15, an = —5, ai2 = —3, G121 = 2, 022 = 1. Bottom: the elliptic case for ST2, with an = 7 and 
all other parameters as in the saddle case. Note, that the homoclinic bifurcation does not terminate 
at ST2- Instead it terminates on saddle-node line SN. 



in the invariant cci-direction. The coordinates and the seven parameters are related 
by three continuous symmetries: 
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(x 1 ,au,a 2 i,e) m> (Xxi, ja n , \a 2 i, Ae) 
(x2, 012,022) t-t (^2,^12,^022) (1.2) 
(61,62, an, a,2i,a>i2, o>22,e,t) 1 ^ (k6i, k6 2 , Kan, ko 2 i, tai2, ko 2 2, ^e, 

for any A,/j,k ^ 0. In the following, we will use 62 and e as bifurcation parameters, 
fixing an and ai2 to distinguish the topological^/ different bifurcation diagrams. The 
remaining parameters are fixed to b\ = 15, 021 = 2 and 022 = 1- In Figure fTTTI two 
different bifurcation diagrams are shown. In both diagrams the single zero eigenvalue 
interaction (labelled STi) and the double zero eigenvalue interaction (labelled ST2) 
occur. 

The system (|1.1[) has at most four equilibria depending on the parameters. Two 
equilibria are sitting on the xi-axis which is invariant. In Figurc fTTTl we see two saddle- 
node bifurcations (labelled SN). The first saddle-node bifurcation, which is a vertical 
line in both figures, is a collision between equilibria that lie on the xi-axis. The other 
saddle-node bifurcation curve involves the other two equilibria. We also have trans- 
critical bifurcation curve (labelled TC) which occurs when an equilibrium crosses 
the Si-axis. Additional codimension-one bifurcations also appear such as Hopf bi- 
furcation curve (labelled HB), heteroclinic connections (labelled Het) and homoclinic 
loops (labelled Horn) . Continuing further those codimension-one bifurcation curves we 
obtain codimension-two bifurcation points such as a Bogdanov-Takens (BT), saddle- 
node/heteroclinic bifurcation (SNHet) and homoclinic/heteroclinic bifurcation(To). 
We will focus on the description of the dynamics around STi and ST2. The latter 
interaction point organises part of the bifurcation diagram. 

The equilibria on the invariant axis are called the predator-free equilibria. De- 
pending on the parameters, one of the following situations is realised: two predator- 
free equilibria, of which at most one stable, a unique predator-free equilibrium of the 
saddle-node type or the absence of a predator-free equilibrium. The coexistence of 
predator-free equilibria is a consequence of the introduction of constant rate harvest- 
ing or migration, which breaks the invariance of the a^-axis so that the origin, which 
represents the total extinction equilibrium, is shifted along the xi-axis. In addition, 
there are equilibria at which both species survive and these can coexist with predator- 
free equilibria. A second consequence of the introduction of constant rate harvesting 
or migration is the existence of limit cycles, proven to be absent in the original Lotka- 
Volterra model |SJ. The limit cycle is either the sole attractor in the first quadrant 
or it is the boundary of the domain of attraction of an equilibrium with coexisting 
species. The limit cycle can be destroyed in two different ways: either in a saddle-node 
homoclinic bifurcation or in a heteroclinic loop. In the former case we see a time series 
that shows short excursions from a predator-free equilibrium and in the latter case we 
see the population densities alternating between two predator-free equilibria, inter- 
spersed with excursions into the region of coexistence. Obviously, the inclusion of the 
migration or harvesting parameter significantly changes the Lotka-Volterra dynamics. 

2. A single zero eigenvalue. Figure [2~T1 shows the dynamics around the single 
zero eigenvalue interaction STi. Three equilibria are involved in this interaction, 
one of which lies on the invariant axis while the others are created in a saddle-node 
bifurcation. 

2.1. The minimal model. A simple model for the qualitative behaviour shown 
in Figure 12.11 is given by 



x = ax + bx 2 + ex 3 , 



(2.1) 
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where e = ±1. Note, that we can restrict our analysis to the case e = 1, which is related 
to the case e = — 1 through the transformation (x, a, b, t, e) — > (— x, — a, b, — £, — e). 
Also, note that this is the normal form of the transcritical bifurcation extended with 
a third-order term. This model with e = 1 has three equilibria, denoted by 

xq = 0, with a zero eigenvalue iff a = 0, 

X\ = — — + —\/b 2 — 4a, with a zero eigenvalue iff a = and b > or a = — , 

x 2 = — — — — \Jb 2 — 4a, with a zero eigenvalue iff a = and fe < or a = — . 

If we set f(x,a) — ax + bx 1 + x 3 it is straightforward to check the nondegeneracy 
conditions of the saddle-node bifurcation at a = b 2 /4: 

-§£(».*/«>--§ »° d gfe.»V4) = g(^/4) = - 6 , 

(2.2) 

and those of the transcritical bifurcation at a = 0: 

|^(0,0) = 0, A^(0,0) = 1 and 0(0,0) = 26, (2.3) 

from which we can conclude that, in the plane of parameters a and b, a saddle-node 
bifurcation takes place along the line a = 6 2 /4 and a transcritical bifurcation takes 
place along the line a — 0. The only point at which these bifurcations are degenerate 
is the origin, at which only one equilibrium exists. 

2.2. Relation to the cusp normal form. The simple translation 

z = x +l ( 2 - 4 ) 

transforms the minimal model (|2.f j) into the standard unfolding of the cusp bifurcation 

z = fi + vz + z 3 (2.5) 

with unfolding parameters /i and v which are functions of the model parameters a 
and b: 
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= 4,(a,b) = l 3 a _\J ) (2.6) 
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Fig. 2.2. Illustration of the saddle-node-transcritical bifurcation as a nonversal unfolding of 
the cusp bifurcation, a: Unfolding i2.5V with the cusp bifurcation at the origin. The dotted curves 
denote isolines of positive and negative b in model i2. IV . h: Likewise with isolines of positive and 
negative a. c; Bifurcation diagram of model i2. IV with the corresponding lines of constant a, b. The 
transversal intersections B, C and D correspond to saddle-node bifurcations whereas the tangency A 
corresponds to a transcritical bifurcation. 



Thus, we can consider the minimal model of this saddle-node-transcritical interaction 
as an unfolding of the cusp normal form. This unfolding is, however, nonversal because 
the map <fi is non-invertible along part of the bifurcation set. The bifurcation set of 
the cusp unfolding has one component, the well known A-shaped curve of saddle-node 
bifurcations given by 

1^ + 1^=0 (2.7) 

The preimage of this set under <j) consists of two components, given by 

a = —b 2 , at which det(Ddi) = — b 2 and 
4 ' v V! 12 

a = 0, at which det(D0) = 0. (2.8) 

With the exception of the codimension two point at the origin, the map <fi is invertible 
along the first component, corresponding to a saddle-node bifurcation. In contrast, 
the Jacobian of the map has rank one along the second component, which explains 
why this curve corresponds to the more degenerate transcritical bifurcation. In Figure 
12.21 the two bifurcation sets are shown along with lines of constant a and b. A line 
along which b is constant is mapped onto a straight line in the plane of parameters [i 
and v. This line intersects the A-shaped bifurcation set twice, once transversely and 
once in a tangency. A line along a is constant is mapped onto a curve which either 
has no intersection with the bifurcation set (a < 0), has two transversal intersections 
(a > 0) or coincides with the bifurcation set (a = 0). 

2.3. Equivalence to the MLV model. The minimal model (|2.ip is equivalent 
to the reduction of the MLV model (jl.ll) to the one-dimensional centre manifold at 
the saddle-node-transcritical interaction STi. This codimension-two point is located 
at 



h* - 
o 2 — 



6ia 22 a 2 i „ b* 2 



x,=- 2 



D 2 ' a 2 i 

bla 22 D 1 



D 2 



x ;=0, (2.9) 
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where we have defined 

Di = a\\ai2 - 012021, D 2 = 2an022 - 012021. 
After an initial transformation given by 

022 



(2.10) 



(121 



-Z2 + Z 3 , 



e 

%2 = X*2 + Z2 , b = b* + Z4 

the MLV model can be written as the extended system 



* . M12O21 
e H - z 3 , 



D 2 



(2.11) 



b 1 a 12 a 2 i 2 D, 
zi = zi + anz 1 

L>2 2 1 



3 O22 9 2 

2:12-2 H — — -^i 2 2 + 011Z3 + 2a\\Z\Z 3 



-D3 

021 



-Z 2 23 



«22 
021 



-^2^4, 



i 2 = a 2 lZiZ 2 + 02122^3 + -22^4, 
Z3 = 0, 

Zi = 0, 
where we have defined 

D3 = 2aiia 22 - 012021 - 022«2i- 



(2.12) 



(2.13) 



This system has a three-dimensional centre manifold which can be represented locally 
as the graph of a function z\ = ip(z2, z 3 , Z4,). The Taylor expansion of this function is 
found to be 



^2,23,2:4) 



Do 



( CL22D1 2 2 ^3 a 22 

— 5 — z 2 + anz 3 z 2 z 3 H Z2Z4 + h. o. t. (2.14) 

oia 2 iai2 V a 2i a 2i a 2 i 



where "h.o.t." stands for higher order terms. Thus, we find for the dynamics in the 
centre manifold that 



l , s . D 2 z 2 / a 22 D 1 2 2 D 3 

Z2 = {ZA + a2lZ 3 )Z2 + - — n — z 2 + anz 3 Z 2 Z 3 

0lOi2 V a 21 "21 

Now if we scale the dependent variable as 



Q22 
021 



-z 2 z 4 



(122D1D2 



6101202! 



z%, 



-h. o. t. (2.15) 



(2.16) 



we find equation (12.11) with 



/ a 2 2DiD 2 \ 
= sign — , 

V M12 / 

. O11D2 2 , 
O = 02lZ 3 + - Z 3 + Z 4 , 



6lOl2 



b = 



eo 2 i 



a 22 DiD2 



6lOl20 21 



( ° 3 \ 
(z 4 z 3 ). 

022 



(2.17) 



The latter relations define a map between the parameters z 3 and 24 and the parameters 
o and b which is smooth and invertible on an open neighbourhood of the codimension- 

two point (2:3,2:4) = ( j, ia ^ 2 2 n21 ( e — e*) 5 6— b*^j = (0,0). In this computation, we have 

assumed that D± and D2 are not equal to zero to avoid higher order degeneracies. 
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3. A double zero eigenvalue. In Figures [3~T1 and [3721 the bifurcations around 
the saddle-node-transcritical interactions with two zero eigenvalues are shown. Again, 
three equilibria are involved but in this case cycles and connecting orbits are generated. 

3.1. The minimal model. A simple model of this interaction is given by 
x = fx(x,y,a,b) = y, 

V = f2{x 1 y, a, b) = ax + k\by + bx 2 + k2xy + x 2 y + ex 3 + k^x 4 , (3-1) 

where ^1,^2,^3 ^ 0, ^ ^ 2^/2 and e = ±1. For a = b = this model coincides 
with the normal form of a degenerate Bogdanov-Takens bifurcation [6j [TT] . This is a 
codimension three singularity and its versal unfolding has three parameters. The local 
bifurcations present in the versal unfolding are saddle-node and Hopf bifurcations of 
equilibria and saddle-node bifurcations of cycles. Just like in the case of the single 
zero saddle-node-transcritical interaction of Sec. [5J our model is a nonversal unfolding 
which induces transcritical bifurcations. 
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In the unfolding of the DBT singularity there are also heteroclinic connections to 
equilibria. Generically, these are one-way connections, in contrast to the heteroclinic 
loops we observe in the MLV model (see Fig. 13. ip . The cause of this structural 
difference is a special property of the MLV model which is not automatically preserved 
in the minimal model. Up to two equilibria of the MLV model are forced to lie on the 
invariant axis. If both are of saddle type, a structurally stable heteroclinic connection 
exists in the model. We can keep this structure in the minimal model if we impose 
some conditions on the coefficients. 

Proposition 1. Under the conditions 

2ek\ - kxk 2 -1 = 

3fcifc 3 -1 = (3.2) 

the manifold given by 

y = g(x, a, b) = ak\ + bk\x + ek\x 2 + —x 3 (3-3) 

is invariant in system i3.1\) . Moreover, all equilibria except the origin lie on this 
manifold. 

Proof: We have f 2 (x, y, a, 6) = xP(x, a, b) + yQ{x, 6), where P and Q are polyno- 
mials in x of order 3 and 2, respectively, and g(x,a,b) = kiP(x,a,b). The manifold 
is invariant if 

g'(x,a,b)fi(x,y,a,b) = f 2 (x,g(x),a,b) 
k 2 P'(x, a,b) = x + kiQ(x, b) 

and this equation holds identically if and only if conditions (|3.2j) are satisfied. The 
observation about the equilibria follows directly from the fact that the equilibria are 
given by y = and f 2 (x, 0, a, b) = xg(x, a, b)/k\ = 0. □ 
In Sec. 13.51 we will show that the MLV model can brought to the form of our 
minimal model by a normal-form transformation. If we compute the corresponding 
transformation in parameter space, we find that conditions (|3.2|) are identically sat- 
isfied. The bifurcation diagrams which arise in the minimal model without these 
conditions are numerous and rich. A complete description falls outside the scope of 
the present paper and will be presented elsewhere. In the following, we will assume 
that conditions (|3.2p hold. 

3.1.1. Basic bifurcation structure. Note, that the model is invariant under 
the reflection (x, y, a, b, k\, k%) — > (— x, — y, a, — b, — k%, — &2, — £3). As a consequence 
we can restrict the description of the bifurcation diagrams to the case k 3 > 0. 
The equilibrium solutions are: 

(xq, 0) = (0, 0), with a zero eigenvalue along the line TC given by a TC = 0, 

(;ci,0) and (x2,0), where Xi t 2 are the roots of f 2 (x,0,a,b) — which coincide with 
xq in the limit of a, 6 — > and have a zero eigenvalue along the line SN given by 



(x3, 0), where X3 is the root of f 2 (x, 0, a, b) = which tends to —e/k^ in the limit of 
a, b-> 0. 
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The latter equilibrium does not play a role in the unfolding of the saddle-node- 
transcritical interaction. 

In order to check the nondegeneracy conditions along the lines TC and SN we 
computed the parameter-dependent centre manifold reductions. Along TC the dy- 
namics on the centre manifold of the origin is given up to third order in x and a 

by 

k\bx — —bx 2 — ax 

which is, up to a scaling, the normal form of the transcritical bifurcation. Along SN, 
the bifurcating equilibrium is located at (x SN ,0) = {eki[yjl — 3bk 3 — 1],0) and the 
dynamics on its centre manifold is given up to third order in x and a by 

x = ki(a — a SN ) + efciyl — 3k 3 b(x — x SN ) 2 

which is, up to a scaling, the normal form of the saddle-node bifurcation. Thus, we 
conclude that equilibria (xi,0) and (x2,0) coalesce in a nondegenerate saddle- node 
bifurcation along SN and either of them cross equilibrium (xq , 0) in a nondegenerate 
transcritical bifurcation along TC. 

In addition to the transcritical bifurcation, the equilibrium (xq , 0) undergoes a 
Hopf bifurcation along the line HB given by b = 0, a < 0. For b — and a > 
this equilibrium is a neutral saddle. Along the HB the Lyapunov coefficient is strictly 
positive, so the bifurcation is nondegenerate away from the codimension two point. 

At the codimension two point a = b — 0, the minimal model ()3.1[) coincides with 
the normal form of a DBT bifurcation for which the topological phase portraits have 
been categorised as follows (see Dumortier et al. [B]): 

• for e = 1 the origin is a topological saddle, 

• for e = — 1 the origin is a topological focus if k\ — 8 < , 

• for e = — 1 the origin is a topological elliptic point if fcf — 8 > 0. 

We will only consider the saddle case and the elliptic case, because the conditions 
(|3.2p imply that k% > 8. Geometrically, this restriction makes sense as the invariant 
manifold given by (|3.3j) passes through the origin for a = b = so it cannot be a 
topological focus. 

3.2. Unfoldings of the saddle case. In Fig. 13.31 the unfolding of the saddle- 
node transcritical interaction is shown for the saddle case. Note, that in the left 
half plane a structurally stable heteroclinic connection between two saddle points 
exists, as explained above. Starting from region 1 and going around in a clockwise 
direction, we first see a Hopf bifurcation of the origin. The cycle grows and becomes 
a heteroclinic cycle on the line Het. After that, one of the saddle points crosses 
the origin in a transcritical bifurcation and becomes a sink. It subsequently collides 
with the remaining saddle in a saddle-node bifurcation. When crossing the line of 
saddle-nodes again, a saddle and a source are created on the other side of the origin. 

3.3. Unfoldings of the elliptic case. In Fig. 13.41 the unfolding of the saddle- 
node transcritical interaction is shown for the elliptic case. Starting from region 1 and 
going around in a clockwise direction, we first see a Hopf bifurcation of the origin. 
The cycle grows and becomes a homoclinic loop to the saddle-node which exists along 
SNo. On TC, the origin changes from a sink to a saddle. The second time we cross 
TC, both equilibria have moved to the left of the origin, after which they collide at 
SN. 
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Fig. 3.3. Bifurcation diagram of the saddle-case of the system i3. It under conditions i3.2t . 
The curve C is a clockwise walk around the codimension two point. Its correspondence to a path in 
the parameter space of the DBT unfolding is shown in Fig. \3.5\ 




Fig. 3.4. Bifurcation diagram of the elliptic-case of the system f3. lt l with conditions $3. 2( 1. 



3.4. Relation to the DBT normal form. We have discussed that at the point 
a = b = 0, the system (|3.1[) becomes the normal form of the degenerate Bogdanov- 
Takens bifurcation. In an open neighbourhood of this point we can define a transfer- 
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mation which relates the two. It is given by 

zi=x+ l b -w 2 xb -2k b2 > 

3fc 2 

a = a- |& 2 , 

l = b-^\ (3.4) 

For the new variables we find the standard unfolding of the DBT bifurcation, trun- 
cated up to terms of order three: 

z\ = z 2 , 

z 2 = Ml + M2^i + vz 2 + k 2 z x z 2 + z\z 2 + tz\, (3.5) 



where 



Mi = ~|( a + ' 6 ' 
M2 = a, 



v 



(h - |fe 2 ) b. (3.6) 



Equations (|3.4I) and (|3 . 6[) define a map $ from the two-dimensional space of parame- 
ters a and b to the three-dimensional space of parameters /Lti, /i 2 and Let us denote 
by S the embedding of some open neighbourhood of the saddle-node-transcritical 
point under <!>. This surface is given by 

[3fci - efe 2 ] 3 Mi + e[3fci - tk 2 ] 2 ^ 2 v + u 3 = 

Also, let SN denote the saddle-node surface of the DBT normal form, given by 

27 'n\ + 4ef4 = 

The surfaces 5* and SN intersect transversely along the curve T SN given by 

Mi = -|[3fc! -efc 2 ]~V // 2 = -|[3fci-efc 2 ]- 2 i/ 2 (3.7) 

and have a tangency along the curve T TC given by 

fix = 2[3fci - efc 2 ]~V /i 2 = -3[3fci - efc 2 ] _ V. (3.8) 

The curves r SN and T TC are the image under $ of the saddle-node and transcritical 
bifurcation lines of the minimal model, respectively. Diagram 13.51 shown the embed- 
ding surface, along with the local bifurcations, for the saddle case. In addition to 
the surface SN of saddle-node bifurcations there is a surface of Hopf bifurcations, 
labelled HB. The label NS denotes a neutral saddle which is not a bifurcation. The 
codimension-two Bogdanov-Takens bifurcation (labelled BT) is a curve along the in- 
tersection of the saddle-node and the Hopf/neutral saddle surfaces. The degenerate 
Bogdanov-Takens bifurcation (labelled BT3) is the origin of this parameter space. 
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Fig. 3.5. Top: schematic partial bifurcation diagram of the DBT singularity 13. 51) for the 
saddle case along with the embedding surface \3.4^ - The parameter values are e = 1, k\ = fea = 1. 
The bifurcation diagram is described in detail in \lh Bottom: the saddle-node surface (transparent) 
and the H op f surface (opaque). The two intersect along the line of Bogdanov-Takens bifurcations 
(solid line). The dotted line is the image of a circle around the origin of the parameter space of the 
minimal model f3.1\) as shown in Fig. \3.3\ The numbers correspond to the phase portraits. The 
solid circles denote saddle-node, transcritical and Hopf bifurcations. The solid square denotes the 
heteroclinic bifurcation. 



We have not drawn surfaces of global bifurcations in Fig. 13.31 No explicit ex- 
pressions of these surfaces are known, but their topology is partly proven and partly 
conjectured by Dumortier et al [6] . The structurally stable heteroclinic connection is 
broken by transformation Q3.4p . This is not a consequence of the truncation to third 
order. In order to compare the unfolding of the saddle-node-transcritical bifurcation 
to that of the degenerate Bogdanov-Takens bifurcation, we have to assume that the 
embedding surface S coincides with a surface of heteroclinic connections if two saddle 
points exist on the invariant manifold given by (|3.3[) . An inspection of the unfoldings 
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of the saddle and elliptic cases in reference [5] shows that the bifurcations diagrams 
presented in Sees. 13.21 and 13.31 are the only possible unfoldings. 



3.5. Equivalence with the MLV model. The saddlc-node-transcritical bi- 
furcation with double-zero eigenvalues occurs in the MLV model when 

xi = x l - 



2an 4an 

x 2 =x* 2 = 0, b 2 = b* 2 = (3.9) 

2an 

We introduce u\ = xi — x*, u 2 = xi — x%, pi = e — e* and p 2 = b 2 — b 2 , thus we have 

ui = g 1 (u 1 ,u 1 ,p 1 ,p 2 ) = ju 2 + anu\ + a^u^+Pi, 

u 2 = g 2 {ui,u\,p\,p 2 ) = a 21 u\u 2 + a 22 u\ +p 2 u 2 , (3.10) 

where 7 = — b\a\ 2 / (2an). Now consider a transformation given by 

ui = m - ^-p\ + —p 2 + 4>i(ui,u 2 ,pi,p 2 ) q\=P2 - n Ds Pi + ipi(pi,p 2 ) 
a 2 ij a 21 27011 

v 2 = ju 2 +pi + c/) 2 {ui,u 2 ,pi,p 2 ), q 2 = pi + 4> 2 (pi,p 2 ) 

T = t[l + 6(u 1 ,u 2 , Pl ,p 2 )} (3.11) 

where 4>i^ 2 and ^1,2 are polynomials in all their variables with zero linear and constant 
parts and 9 is a polynomial with zero constant part. Clearly, this transformation is 
smooth and invertible on an open neighbourhood of the codimension-two point. The 
equation for v\ can be normalised by choosing the coefficients of <p 2 so that, up to 
fourth order 



4>2 = [1 + 0] 1 {gi + gid Ul (/)i +g 2 d U2 (/)i) - ju 2 



Pi 



It is a straightforward if tedious exercise to choose the coefficients of </>i, -01,2 and 
9 to normalise the equation for v 2 . From the theory of the DBT singularity we 
know that all terms of the form v^v™ can be removed, except when m = 0, 1, and 
moreover the term vfv 2 can be removed by hypernormalization. The elimination of 
terms involving the parameters requires solving about fifty linear equations with up 
to a few thousand terms, which is best done using a computer algebra system. As the 
resulting transformation contains as many terms, we omit the details. The result up 
to fourth order is the following ODE: 

dvi 

dv 2 2<2n 2 OL11D3 2 o 

— — = -a 2 iq 2 vi qiv 2 + 2anqiv 1 + D4V1V2 H v x v 2 - a2\a\\v 1 

dl a 21 7a 2 i 

160^1)3 2 2 4a u Z) 3 (2an - a 2 i) 3 anD 3 4 
+ — TTT9i w i o n V lV i 5 u i ( 3 - 12 ) 



The saddle-node-transcritical bifurcation 



15 



where D4 = 2an + a 21 . Finally, we scale the variables and time as 



D3J \ana 2 \\ _ 2aiiD 3 
—^2 v i 9i = — 7. — ,, 

7«21 7«2lVl a ll«21 



^lVl a n«2i| _ Dl 

v 2 = -e o~i v 2 92 = 5-3- Q2 

r = - e 2^ii T (3.13) 

where e = — sign(ana2i). This gives 

«1 = «2 

«2 = <Z2^i + kiqiv 2 + qiv\ + k 2 viV2 + v\v 2 + ev\ + k 3 vj + k±q~lvl + k 5 qivf (3.14) 
where 



, £\/l a ii a 2il , e -°4 , \f\au02i\ 

ki = k 2 = ,, k 3 = 



a>2i \/l aiia 2il 



1 



_ 16a nV /|aiia 2 i| _ 4e(2an - Q2i)\/|Qiia2i| ,„ , 

fc4 ~ 3^ % _ 3Dl 

Note, that fci, fc 2 and £3 identically satisfy conditions (|3.2[) . The terms proportional 
to A4 and £5, however, introduce a splitting of the invariant manifolds of the saddle 
type equilibria. Effectively, these terms unfold the heteroclinic loop of the MLV model 
into two separate heteroclinic connections and a homoclinic bifurcation. A complete 
description of the unfolding of the saddle-node-transcritical interaction in the absence 
of the special structure of the MLV model, reflected by Proposition [TJ is out of the 
scope of the present paper and will be presented elsewhere. For our present purpose 
it suffices to simply neglect the extra terms, in which case Eqs. (|3.14[) is identical to 
the minimal model (|3.1[) . 

The nondegeneracy conditions for the transformation are 

7)°ii> 0121)022) D 3 ,D A , (a 2 2 + ai 2 ), (an - 021), (4on - a 2 i) 7^ 

4. Conclusion. Standard codim-2 bifurcations, such as cusp and Bogdanov- 
Takens bifurcations, have been widely investigated in mathematical models of popula- 
tion dynamics (see, e.g. [T7J 1131 QUI [H] ) ■ In this paper, we investigated a non-standard 
codimension-two bifurcation, namely the interaction of saddle-node and transcritical 
bifurcations, in a Lotka-Volterra model modified to describe harvesting or migration. 
We have shown that the two different interactions, associated to either a single zero 
eigenvalue or a pair of zero eigenvalues, are described by nonversal unfoldings of the 
standard cusp and degenerate Bogdanov-Takens bifurcations. In the latter case so- 
lutions exist which are not allowed in the original Lotka-Volterra system, namely 
periodic, heteroclinic and homoclinic solutions. 

Is is somewhat surprising that a small modification of the Lotka-Volterra model 
so significantly changes its dynamics. In the modified model, we see coexistence of two 
predator-free equilibria and periodic fluctuations of the densities of coexisting species. 
These fluctuations can have arbitrary long periods and model short excursions from 
predator-free states. 
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From a mathematical point of view, work that needs to be done includes the 
analysis of the minimal model (|3.ip in the absence of the special structure which 
imposes the presence of a structurally stable heteroclinic connection. This is work in 
progress and will presented elsewhere. We hope that the description of the interaction 
of the transcritical bifurcation with other local bifurcations will yield new tools to 
analyse models in which such interactions are of codimension two. 
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